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ABSTRACT 



The peculiar velocity function of clusters of galaxies is determined using 
an accurate sample of cluster velocities based on Tully-Fisher distances of Sc 
galaxies (Giovanelli et al 1995b). In contrast with previous results based on 
samples with considerably larger velocity uncertainties, the observed velocity 
function does not exhibit a tail of high velocity clusters . The results indicate 
a low probability of ^ 5% of finding clusters with one-dimensional velocities 
greater than ~ 600 km s~^. The root-mean-square one-dimensional cluster 
velocity is 293±28 km s^^. The observed cluster velocity function is compared 
with expectations from different cosmological models. The absence of a high 
velocity tail in the observed function is most consistent with a low mass-density 
(n ~0.3) CDM model, and is inconsistent at ^ 3a level with Q= 1.0 CDM and 
HDM models. The root-mean-square one-dimensional cluster velocities in these 
models correspond, respectively, to 314, 516, and 632 km s~^ (when convolved 
with the observational uncertainties). Comparison with the observed RMS 
cluster velocity of 293±28 km s^^ further supports the low-density CDM model. 

Subject headings: galaxies:clusters:general — cosmology:observations — 
cosmology:theory — dark matter — large-scale structure of universe 
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1. Introduction 



The motions of clusters of galaxies can place strong constraints on cosmological 
models and on the mass-density of the universe. Bahcall et al (1994a, b), Cen et al (1994), 
Croft & Efstathiou(1994), Lauer & Postman(1994), Gramann et al (1995), and Moscardini 
et al (1995), showed that clusters of galaxies provide a particularly efficient and accurate 
way to trace the peculiar velocity field in the universe. In turn, the peculiar velocity field, 
caused by the gravitational growth of structure, sheds light on the cosmology responsible 
for the formation and evolution of the structure (Dekel 1994, Strauss & Willick 1995). 

Bahcall et al (1994a, b) investigated the probability distribution function of cluster 
peculiar velocities, i.e., the cluster velocity function (CVF), and showed that it provides an 
important tool for distinguishing between different cosmological models. They determined 
the cluster velocity function for several cosmological models using large scale N-body 
simulations. They also determined the observed CVF using the available data and compared 
it with model expectations. However, the large uncertainties of the cluster velocity data 
broadened the CVF and produced an artificial tail of high velocity clusters. These 
uncertainties did not allow an accurate determination of the true underlying cluster velocity 
function, nor an accurate comparison with the cosmological models (since the convolution 
of the model CVFs with large observational uncertainties reduced the differences between 
the various models). Similar results were also obtained by Croft &: Efstathiou(1994) and 
Moscardini et al (1995). Bahcall et al (1994b) concluded that a cluster sample with 
considerably improved velocity accuracy is needed before an accurate cluster velocity 
function, one that is not dominated by velocity errors, can be determined. 

In this paper, we use a new sample (Giovanelli et al 1995a, b) of cluster velocities that 
has considerably higher accuracy and uniformity than previously used samples. The new 
sample is based on well calibrated TtiUy-Fisher distance indicators of Sc galaxies. We use 
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these data to determine the cluster velocity function and compare it with expectations from 
cosmological models. 

2. The Peculiar Velocity Function of Clusters 

2.1. Model Expectations 

The peculiar velocity function of clusters of galaxies represents the probability 
distribution of cluster velocities relative to a comoving cosmic frame. The integrated 
velocity function, P(>v), represents the relative number density of clusters with peculiar 
velocities larger than v (where v is the three-dimensional cluster motion relative to the 
cosmic frame). The differential velocity function, P(v), represents the relative number 
density of clusters with peculiar velocities in the range v±dv, per unit dv, as a function of 
V. The cluster velocity functions P(v) and P(>v) were determined for four cosmological 
models by Bahcall et al (1994b) using large scale N-body simulations. We use these results 
below. 

The cosmological models investigated and their parameters are summarized in Table 
1. These parameters include the matter density, fl; the cosmological constant contribution, 
^Aithe Hubble constant (in units of Hg = lOO/i km s~^ Mpc~^); and the normalization of 
the mass fluctuations on a 8h^^ Mpc scale, cjg. The models are normalized to the large-scale 
microwave background anisotropy measured by COBE (Smoot et al 1992). (The HDM 
model normalization is ~ 20% higher than the il= 1 CDM on large scales). We next 
describe briefly the simulations that are used to represent the cosmological models. 

A large-scale particle-mesh code with box size of 800 Mpc is used to simulate the 
evolution of the dark matter in the models. A large simulation box is needed in order to 
ensure that contributions to velocities from waves larger than the box size are small, and 
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to minimize uncertainties due to fluctuations in the small number of large waves. The 
simulation box contains 500^ cells and 250^ = 10^'^ dark matter particles. The spatial 
resolution is 1.6 h^^ Mpc. (A higher resolution [0.8 h^^ Mpc] smaller box [400 h^^ Mpc] was 
also studied for comparison). For more details of the simulations see Bahcall et al (1994b). 

Clusters were selected in each simulation using an adaptive linkage algorithm. The 
cluster mass thresholds correspond to the observed number density of typical rich clusters 
as well as of groups. A total of ~3000 rich clusters and ~5 x 10^ groups were obtained 
in each of the simulated models. The three-dimensional and one-dimensional peculiar 
velocity of each cluster or group, relative to the comoving cosmic frame, was obtained 
from the simulation and used to determine the velocity function of groups and clusters. 
The simulation results are consistent with expectations from linear theory (Bahcall et al , 
1994a, b). The clusters selected for comparison with the present data correspond to the 
group selection threshold, which represents the best match to the threshold of the observed 
groups and clusters in the current sample. The results, however, are insensitive to the exact 
richness threshold of the clusters (Bahcall et al 1994b). 

The cluster velocity functions of the four models are presented in Figures 3-4 and 9-11 
of Bahcall et al (1994b); these functions represent the "exact" CVFs (in vsd and yio), 
unconvolved with any observational uncertainties. The results illustrate that the diff'erences 
among the four models are most apparent at the high velocity end, where the low-density 
models predict considerably smaller peculiar velocities than the f]= 1 models. For example, 
while the Q,= 0.3 CDM and FBI models yield ^ 5% of clusters with velocities vm > 500 
km s^^ and > 800 km s^^ respectively, the = 1 CDM and HDM models exhibit ~5% 
of clusters with high velocities of Vid > 1000 km s~^ and > 1300 km s~^, respectively. 
Similarly, the root-mean-square peculiar velocity of clusters diff'ers significantly among 
the models. The n=0.3 CDM model yields the lowest RMS velocity, < v^^ >^/^ ~ 268 
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km s^^, while the Q=l models yield the highest velocities, < v^jj >^/^ ~ 500-600 km s^^ 
(unconvolved with observational uncertainties). The results are summarized in Table 3 of 
Bahcall et al (1994b). The sensitivity of the CVF to the cosmology makes it a powerful 
tool in constraining cosmological models. We use this tool below. 

2.2. Observations 

The first determination of the cluster velocity function was made by Bahcall 
et al (1994b) who used observations of cluster velocities based on Tully-Fisher (TF) and 
D„ — (7 distance indicators (with data from Aaronson et al 1986, Faber et al 1989, Mould 
et al 1991, 1993, and Mathewson, Ford & Buchhorn 1992). They found a velocity function 
that exhibits a large tail of high velocity clusters up to yib ~ 2000 km s^^ (Figs. 10-11 
of Bahcall et al 1994b). However, the observational uncertainties of the cluster velocities 
were very large, reaching ~900 km s~^. The authors showed that when the model cluster 
velocities are convolved with the large observational velocity uncertainties, an artificial high 
velocity tail, not present in the original model CVF, is produced. Even this artificial high 
velocity tail was in general not as large as suggested by the data (especially the D„ — a 
data). Bahcall et al suggested that the high velocity tail of the CVF was an artifact of large 
velocity uncertainties. Differences between the observations and model expectations could 
arise from imderestimated velocity errors. The authors emphasized the need for a cluster 
sample with higher velocity accuracy in order to better determine the CVF, especially at 
the critical high velocity end. 

Recently, a uniform and accurate sample of peculiar velocities of clusters was obtained 
by Giovanelli et al (1995a, b). Their cluster velocities have considerably greater accuracy 
than previous studies, mainly due to: (a) access to a homogeneous, all-sky survey, (b) a 
different TF template relation, based on an extensive study of clusters, (c) an internal 
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extinction correction that allows for larger flux corrections and is luminosity dependent. 
While the sample size is small (22 groups and clusters out to cz < 10 000 km s~^), which 
can thus introduce significant statistical uncertainties, the high quality of the velocity 
measurements provides a clear advantage. The cluster velocity uncertainties range from 
~ 50 km s~^ to 340 km s~^, with a mean uncertainty of 160 km s~^. By contrast, the 
previous sample uncertainties ranged from ^ 70 km s^^ to ^900 km s^^, with a mean 
uncertainty of 410 km s^^. We use this sample to determine the CVF and to compare it 
with model expectations. The reduced observational uncertainties increase the accuracy of 
the measured CVF, especially at higher velocities. 

We present in Figs, la-c the CVF determined from the Giovanelli et al (1995a, b) 
sample. (The 1995b sample is slightly larger, with some improvements over the 1995a 
sample; both overlapping samples yield similar results, and both are presented in Fig. la, 
for comparison.) The error-bars correspond to ±y/N statistical uncertainties. The curves 
in Fig. 1 represent the CVFs of the four cosmological models (§2.1) convolved with the 
observational velocity uncertainties (for proper comparison with the data). The small 
velocity uncertainties of this sample have only a minor impact on the true (unconvolved) 
CVF. 

The new data, in contrast with previous samples, do not exhibit a high velocity tail. 
In fact, there are no observed clusters with velocities larger than vm ~600 km s^\ yielding 
P(viD > 600 km s~^)^ 0.05. In contrast, the previous CVF based on data with larger 
velocity uncertainties showed a high velocity tail to vid ^ 2000 km s~^, with P{vid> 
600 km s-i)- 0.4, and P{vid> 1000 km s^^) -0.1 (Bahcall et al 1994b). Similarly, the 
root-mean-square velocity of the current cluster sample is < v^^, >^/^ = 293±28km s~\ 
as compared with 607±64 km s~^ for the previous TF data and 725±50 km s~^ for the 
previous TF + D„cr data (Bahcall et al 1994b). 
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3. Comparison of Models with Observations 



The observed and model velocity functions are compared with each other in Figure 1 
and Table 2. The main difference among the models is apparent: the Q=l models (CDM 
and HDM, convolved with the observational velocity uncertainties) exhibit a large tail of 
high velocity clusters, with ~5% of all clusters at vm > 1100 km s^^ and 1450 km s^^ 
respectively, while fl = 0.3 CDM exhibits the lowest cluster velocities, with ~5% of clusters 
at vir, > 650 km s~^ (for the convolved models), il = 0.3 PBI has intermediate velocities, 
with ~5% of clusters at vm > 900 km s^\ The observed CVF indicates a clear absence of 
high velocity clusters. This is consistent with the D, ~ 0.3 CDM model and inconsistent 
(at ^ 3a) with the 0,= 1 CDM and HDM models. We do not observe any clusters with 

>600 km s^^ in this sample; the observed CVF yields P(vi£) > 600 km s^^)^ 0.05, or 
^ 1 cluster out of a sample of 22 clusters. From the integrated model CVFs (Figs. lb,c) 
we would expect on average to find 1.7, 4, 6 and 8 clusters with vid >600 km s~^ in a 
random sample of 22 clusters for Q = 0.3 CDM, Q = 0.3 PBI, Q = 1 CDM, and 1^ = 1 
HDM respectively (Table 2). The probability that the observed CVF for vm > 600 km s"^ 
is consistent with the various models can be estimated using the binomial distribution 
statistic, yielding significance levels of 48%, 4%, 1%, and < 0.1%, for Q = 0.3 CDM, n = 
0.3 PBI, n = 1 CDM, and ^] = 1 HDM respectively. A formal K-S test of the integrated 
CVF (Fig. 2b) indicates that the data is consistent with the models at significance levels of 
- 90%, 13%, 1% and 0.1% respectively for fl = 0.3 CDM, fl = 0.3 PBI, fl = 1 CDM, and 
n = 1 HDM. A mixed dark matter model, with -70% CDM and -30% HDM, is expected 
to yield results similar to Q — 1 CDM. The results are summarized in Table 2. 

The RMS peculiar velocity of clusters in the present sample (§2.2) is compared 
with model expectations in Table 2. The RMS velocity, as well as the K-S test, and an 
inspection of Figures la-c all suggest that the fl= 0.3 CDM model is consistent with the 
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observed cluster velocities, while the 0,= 1 models are less consistent with the data. (A 
PBI model with fl ^ 0.3 is also acceptable). The present sample allows, for the first time, 
an accurate comparison of cluster velocities with model expectations, not dominated by 
velocity uncertainties. A larger sample, with similarly accurate velocities, is needed in order 
to confirm and refine this conclusion. 

4. Conclusions 

We have determined the peculiar velocity function of clusters of galaxies using a 
small but accurate sample of cluster velocities (Giovanelli et al 1995a, b). The relatively 
accurate velocities enable a reliable determination of the CVF, not dominated by velocity 
uncertainties. The CVF shows no high velocity clusters, P(vi£) >600 km s^^);^ 0.05, in 
contrast with less accurate previous samples that exhibited an artificially large velocity tail 
to VjQ ~ 2000 km s~\ with P{vid >600 km s~^) ~ 0.4. The root-mean-square ID cluster 
velocity is < vf^ >i/2= 293±28 km s'K 

We compare the cluster VF with expectations from several cosmological models. We 
find the data to be most consistent with a low-density {Q ~ 0.3) flat CDM model, marginally 
consistent with a low-density flat PBI model (Q ~ 0.3), and inconsistent at ^ 3a level with 
Q= 1 CDM and HDM models in which a larger high velocity tail is expected. Similarly, 
the RMS cluster ID velocities in the models yield (when convolved with the observational 
uncertainties) 314, 423, 516, and 632 km s^^, respectively, as compared with the observed 
293 ± 28 km s"^ further supporting the ^ 0.3 CDM model. A low-density flat CDM 
model, which best fits other observations, including the mass function and correlation 
function of clusters (Bahcall & Cen 1992), the baryon density in clusters (White et al 1993, 
Lubin et al 1996), the power spectrum and small scale velocities of galaxies (Maddox et al 
1990, Ostriker 1993), is therefore also consistent with the cluster velocity function. 
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5. Figure Captions 



Fig. 1. — Observed versus model cluster velocity functions. The Giovanelli etal (1995a, b) 
IRTF velocity data(§2.2), is compared with model CVFs convolved with the observational 
errors. (Due to the small observational errors, the effect of the convolution is small). Fig la 
represents the differential function, Figs, lb & Ic represent the integrated function on linear 
and log scales. The solid line histogram (with ^/N statistical error bars) represents the 22 
cluster sample (1995b), while the points (plotted only for Fig. la) represent the earlier 16 
cluster sample (1995a). Note the absence of a high- velocity tail in the observed CVF. 
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Fig. la. — 
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Table 1. Model Parameters 



Parameter 



Model 


n 




h 




CDM 


1.0 


0.0 


0.50 


1.05 


CDM 


0.3 


0.7 


0.67 


0.67 


HDM 


1.0 


0.0 


0.50 


0.86 


PBI 


0.3 


0.7 


0.50 


1.02 



''CDM = Cold Dark Matter model; HDM = Hot Dark Matter model; 
PBI = Primeval Baryonic Isocurvature model. 
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Table 2. Cluster Velocities: Models vs. Observations 





(a) 


(b) 


(c) 




vrms 


Velocity Function 


Nci 




(km s^^) 


KS test 


(viD > 600 km s-i) 


Observed 


293 ± 28 




< 1 


CDM n = 0.3 


314 


90% 


1.7 


PBI n = 0.3 


423 


13% 


4 


CDM n = 1.0 


516 


1% 


6 


HDM ^l = 1.0 


632 


< 0.1% 


8 



^Cluster RMS velocities, < v^j-, >^/^. The model velocities are convolved with the 
observational velocity uncertainties. 

''The KS significance levels for the observed vs. model cumulative velocity functions. 
•^The number of observed versus expected clusters with velocities vm > 600km s^^, 

(for a random sample of 22 clusters for each model) . 
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